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Abstract 

We present results for the negative parity low-lying state of the nucleon, Nj (1535 MeV) Sn, from a variational 
analysis method. The analysis is performed in quenched QCD with the FLIC fermion action. The principal focus 
of this paper is to explore the level ordering between the Roper (Pn) and the negative parity ground (Sn) states of 
the nucleon. Evidence of the physical level ordering is observed at light quark masses. A wide variety of smeared- 
smeared correlation functions are used to construct correlation matrices. A comprehensive correlation matrix analysis 



is performed to ensure an accurate isolation of the A^:^ state. 
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1. Introduction 

Lattice QCD is very successful in computing many 
properties of hadrons from first principles. In particular, 
the ground state of the hadron spectrum is a well un- 
derstood problem |1]. However, the excited states still 
prove a significant challenge. One of the long-standing 
puzzles in hadron spectroscopy has been the low mass 
of the first positive parity excitation of the nucleon, 
known as the Roper resonance, A^^ (1440 MeV)Pii, 
compared to the lowest-lying negative parity partner, 
N^ (1535 MeV)Sii. This phenomenon cannot be ob- 
served in constituent or valence quark models where the 
lowest-lying odd parity state naturally occurs below the 
N - J state. Similar difficulties in the level orderings 
also appear for the J^ = |^A*(1600) and i^5:*(1690) 
resonances. 

There has been extensive research focussing on the is- 
sue of the level ordering problem using the lattice QCD 
approach 021 41411 . Oneof the state-of-the-art approaches 
that has been used extensively in hadron spectroscopy 
is the 'variational method' il5lll6ll . which is based on 
a correlation matrix analysis. In the past the isolation 
of the Roper resonance was elusive with this method. 
However, in Refs. fllTlllSll a low-lying Roper state has 
been identified using a correlation matrix construction 
with smeared-smeared correlators. Our work there mo- 
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tivates us to investigate the long-standing level ordering 
problem using the same techniques on the same lattice. 

In contrast to the positive parity ground state of 
the nucleon, Nj , which has a large plateau over 
Euclidean-time, the correlation functions for the neg- 
ative parity ground state, N^ , are short-lived giving 
shorter plateaus at earlier Euclidean times. Therefore, 
the standard analysis to extract the N^ ground state 
from small Euclidean-times may provide a mixture of 
ground and excited states. On the other hand, the vari- 
ational method accounts for the presence of excited 
states in the correlation functions via correlation matri- 
ces. The masses of the energy states are then obtained 
by projecting the correlation matrix to eigenstates [11711 
providing a robust approach for extracting the energy 
states. In addition, considering several bases in con- 
structing different correlation matrices provides a sub- 
stantial verification of the analysis technique, in allow- 
ing the consistency of the energy states over the differ- 
ent basis and the reliability of the extracted eigenstates 
energies to be explored. 

In this paper, we use the same approach as that of 
Ref. |il71 to isolate the negative parity states of the nu- 
cleon. In particular, we focus on the negative parity state 
to explore the level ordering problem. 
Various sweeps of gauge invariant Gaussian smearing 
lll9ll are used to construct a smeared-smeared correla- 
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tion function basis to form correlation matrices. 

This paper is arranged as follows: Section[2]contains 
a brief description of the variational method. Lattice de- 
tails are in Section|3] Results are discussed in Section|4] 
and conclusions are presented in Section|5] 



2. Variational Method 

The two point correlation function matrix for p - Q 
can be written as 



Gp^ = y Tr,p{ri<Qbr,W;ey(0)|n)), 



Z 



AfA^je-"'"', 



(1) 



(2) 



where, Dirac indices are implicit. Here, A" and A" are 
the couplings of interpolators Xi and Xj at the sink and 
source respectively and a enumerates the energy eigen- 
states with mass m„. Fj- = ^(yo ± 1) projects the parity 
of the eigenstates. 

Since the only t dependence comes from the exponen- 
tial term, one can seek a linear superposition of interpo- 
lators, Xju", such that, 



Gij(to + At)u']^e-"'''''Gijito)u'], 



(3) 



for sufficiently large fp and to + At. More detail can be 
found in Refs. la, |2QI |2l|] . Multiplying the above equa- 
tion by [Gy(fo)]"' from the left leads to an eigenvalue 
equation. 



[(G(?o)r'G(fo + AO],7Mj=c"< 



(4) 



where c" = e^"'"'^' is the eigenvalue. Similar to Eq.(|4]i, 
one can also solve the left eigenvalue equation to re- 
cover the v" eigenvector. 



vf [G(fo + Af)(G(fo))-'],7 



c°v°. 



(5) 



The vectors u" and v" diagonalize the correlation matrix 
at time to and fo + Af making the projected correlation 
matrix, 

vfGf/f)«^ - 6"^. (6) 

The parity projected, eigenstate projected correlator, 

Gl ^ vfGf/O" •, (7) 

is then analyzed using standard techniques to obtain the 
masses of different states. 



3. Simulation Details 

Our lattice ensemble is the same as that explored in 
Ref. II17II . It consists of 200 quenched configurations 
with a lattice volume of 16^ x 32. Gauge field con- 
figurations are generated by using the DBW2 gauge 
action [|22[ 1231] and an (9(a)-improved FLIC fermion 
action ll24ll is used to generate quark propagators. This 
action has excellent scaling properties and provides 
near continuum results at finite lattice spacing 112511 . 
The lattice spacing is a = 0.127 fm, as determined by 
the static quark potential, with the scale set using the 
Sommer scale, vq - 0.49 fm 12611 . In the irrelevant 
operators of the fermion action we apply four sweeps 
of stout-link smearing to the gauge links to reduce 
the coupling with the high frequency modes of the 
theory [27] providing 0{a) improvement |25]. We use 
the same method as in Refs. 1.20. .28.1 to determine fixed 
boundary effects, and the effects are significant only af- 
ter time slice 25 in the present analysis. Various sweeps 
of gauge invariant Gaussian smearing |19] (1, 3, 7, 12, 
16, 26, 35, 48 sweeps) corresponding to rms radii, in 
lattice units, of 0.6897, 1.0459, 1.5831, 2.0639, 2.3792, 
3.0284, 3.5237, 4.1868, are appUed at the source 
(f - 4) and at the sink. This is to ensure a large range 
of overlaps of the interpolators with the lower-lying 
states. The analysis is performed on eight different 
quark masses corresponding to pion masses of m„ - 
{0.797, 0.729, 0.641, 0.541, 0.430, 0.380, 0.327, 0.295) 
GeV The error analysis is performed using the jack- 
knife method, with thej^'^/dof obtained via a covariance 
matrix analysis method. Our fitting method is discussed 
extensively in Ref. |20]. 
The nucleon interpolators we consider are, 

Xiix)^ e"''{u'\x) Cys d\x)) u'(x), (8) 

(9) 



X2(x) = e"'"(u'"(x)Cd'(x))ys u'(x). 



We use the Dirac representation of the gamma matri- 
ces in our analysis. 



4. Results 

We consider several 3x3,4x4, 6x6 and 8x8 corre- 
lation matrices. Each matrix is constructed with differ- 
ent sets of correlation functions, each set element corre- 
sponding to a different numbers of sweeps of gauge in- 
variant Gaussian smearing at the source and sink of the 
XiX\7 X2X2 ^^^ X\X2 correlators |18]. This provides a 
large basis of operators with a variety of overlaps among 
energy states. 
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Figure 1: (colour online). Masses of the nucleon, N^ states, from the projected correlation functions (left) and from the eigenvalues (right), for 
the 3x3 (top row) and 4x4 (bottom row) correlation matrices of ;j'i^i, for the pion mass of 797 MeV. The 3x3 results correspond to the P' 
combination (7,16,26 sweeps) of 3 X 3 matrices, while 4x4 results correspond to the 2 combination (3,12,26,35 sweeps) of 4 x 4 matrices. 
Each set of N^ ground (g.s) and excited (e.s) states masses correspond to the diagonalization of the congelation matrix for each set of variational 
parameters fstaii (shown in major tick marks) and At (shown in minor tick marks). Here, (start - 'o as shown in Eqs. J4) and (5). 



We consider five smearing combinations (bases) 
{1=(7,16,26),2=(7,16,35),3=(12,16,26),4=(12,26,35), 
5 =(16,26,35)} for 3x3 correlation function matrices and 
four combinations {1=(1, 12,26,48), 2=(3,12,26,35), 
3^(3,12,26,48), 4=(12,16,26,35)) for 4 x 4 matrices, of 
X\Xi correlation functions. In the latter case these four 
combinations are found optimal for the reliable extrac- 
tion of the low-lying energy states shown in Ref . 11711 . 
Including the;(f2 interpolator, which vanishes in the non- 
relativistic limit 11291 boll . in correlation matrix anal- 
ysis provides extra challenges. Nonetheless, we con- 
sider this interpolator for the reliable extraction of the 
negative parity ground state mass. The same bases, as 
discussed above for the ;t'i;t'i analysis, are also con- 
sidered for the 3x3 and 4x4 correlation matri- 
ces of xzXi correlation functions. We also consider 
four smearing combinations {1=(3,12,26), 2=(3, 16,48), 
3=(7,16,35), 4=(12,16,26)) of 6 x 6 and four combina- 



tions {1=(3,12,26,48), 2^(7,12,26,35), 3^(7,16,26,35), 
4=(7,16,35,48)) for 8 x 8 matrices of ;ifi;if2 correlation 
functions. 

In Fig.[ri masses from the projected correlation func- 
tions and from the eigenvalues are shown for the 3x3 
and 4x4 correlation matrices. We refer to the low- 
est lying state as the ground state in the negative par- 
ity sector As in Refs. 1171 120(1 . masses from the pro- 



jected correlation functions for the low-lying states are 
very consistent over the variational parameters of fstart 
and Af, in particular, the negative parity ground state is 
robust. However, a deteriorating signal to noise is ev- 
ident for larger fstart and Af values, particularly for the 
excited states. In contrast, the mass from the eigenvalue 
analysis shows significant dependence on the variational 
parameters. Therefore, exposing a mass from the pro- 
jected correlation functions is again proved to be more 
reliable than from the eigenvalues 01711 . 
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Figure 2: (color online). Masses of the nucleon, N^ - states, from projected correlation functions as shown in Eq. for the pion mass of 797 
MeV (left) and 380 MeV (right). Numbers in the horizontal scale coiTespond to each combination of smeared 3x3 correlation matrices of ^1^1 . 
For instance, 1 and 2 correspond to the combinations of (7,16,26) and (7,16,35) sweeps, respectively and so on, as discussed in the text following 
Eq. E). 
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Figure 3: (color online). As in Fig.|2] but for the 4x4 correlation matrices of ^1^1 . 



From a series of fstan and Af, a single mass is selected 
for one set of fstait and At by the selection criteria dis- 
cussed in Ref. 1120(1 . where we prefer larger value of 
f start + Af |2ll]. In cases where a larger fstait + Af provides 
a poor signal-to-noise ratio, for example (fstart, Af)=(7,3) 
(top left graph of Fig. [T]i, we prefer a little lower 
fstart + Af value, for example (fstart, Af)=(7,2), and we fol- 
low this procedure for each quark mass, as discussed in 
Ref. fH. 



In Fig. 121 masses extracted from all the combinations 
of 3 X 3 matrices (from 1'*' to 5*) are shown for the 
pion mass of 797 and 380 MeV. Here the negative par- 
ity ground and the first excited states are very consistent 
for all the 3x3 bases. The second excited state shows 
some smearing dependence [20] which is also observed 
in Ref. lll7ll as this highest excited state must accommo- 



date all remaining spectral strength. 

In Fig. [3] masses extracted from all the combinations 
of 4 X 4 matrices of x\X\ are shown. In Ref. lll7r . the 
three low-lying states were consistent over these four 
bases for the positive parity excited states and similarly 
consistent results are also observed in this negative par- 
ity correlation matrix analysis, again illustrating the ro- 
bustness of the method 1171 . A smearing dependency 
for the highest excited state is also noticed in this 4x4 
analysis, which is to be expected. 

It is noted that although the 3x3 correlation matrix 
analyses are successful for all the eight quark masses, 
the 4x4 correlation matrix analyses are successful only 
for the six heavier quarks. In this case, the signal- 
to-noise ratio for the lighter quarks deteriorates more 
rapidly than the heavier quarks. As a result the vari- 
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Figure 4: (color online). Masses of the nucleon, Nj - states, from 
projected correlation functions as shown in Eq. Q for the pion mass 
of 797 MeV. Numbers in the horizontal scale correspond to each com- 
bination of smeared 3x3 congelation matrices ofxiXl- For instance, 
1 and 2 correspond to the combinations of (7,16,26) and (7,16,35) 
sweeps, respectively and so on. 



Figure 6: (color online). Masses of the nucleon, N^ - states, from 
projected correlation functions as shown in Eq. for the pion mass 
of 797 MeV. Numbers in the horizontal scale correspond to each com- 
bination of sineared 6x6 coirelation inatrices of^i^2- For instance, 
1 and 2 con'espond to the combinations of (3,12,26) and (3,16,48) 
sweeps, respectively and so on. 
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Figure 5: (color online). Masses of the nucleon, N\ - states, from 
projected correlation functions as shown in Eq. (7) for the pion mass 
of 797 MeV. Numbers in the horizontal scale correspond to each com- 
bination of smeared 4x4 correlation matrices ofxiXl ■ For instance, 1 
and 2 correspond to the combinations of (1,12,26,48) and (3,12,26,35) 
sweeps, respectively and so on. 



ational analysis for 4x4 correlation matrices for the 
lighter two quarks is unsuccessful. 

We now discuss the 3x3 and 4x4 correlation matri- 
ces of;t'2A'2- As the;(f2 interpolator is very challenging to 
perform a simulation with and the signal-to-noise ratio 
of the ;t'2^2 correlator deteriorates rapidly, the diagonal- 
ization of 3 x 3 (Fig.|4| and 4x4 (Fig.|5]l matrices are 
only successful for the four heavier quark masses. The 
energy states revealed by the X2 spin-flavour analysis 
are also very consistent as with the XiX\ analysis. In 
particular, the N\ ground state is again robust. 
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Figure 7: (color online). As in Fig. |6] but for the 8x8 correlation 
matrices of^i^2- 



We now consider the 6x6 and 8x8 analyses involving 
both;^fl and;^f2, where linear combinations of ;^i and;t'2 
are utilized to isolate eigenstates. 

The lowest two energy states extracted with the 6x6 
basis are nearly degenerate as illustrated in Fig. |6] An 
explanation for these two nearby degenerate states is 
provided by the quark model. In the quark model based 
on SU(6) symmetry, four decuplet and two octet states 
contribute to provide the 56 representation of SU(6). 
The nucleon and Roper belongs to the ground and ex- 
cited 56 representation of SU(6) symmetry respectively. 
On the other hand, the odd parity ground (1535 MeV) 
and (1650 MeV) states belong to the negative parity, 
L - 1, 70 plet of SU(6). As three spin-| quarks may 
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Figure 8: (color online). Eigenvector values, v" and u", as shown in Eq. (TJ, for the first (lowest) state (top left), second (top right), third (bottom 
left) and fourth (bottom right) energy states. The figure corresponds to the pion mass of 797 MeV and for the 2'"' combination (7, 12,26,35 sweeps) 
of 8 X 8 correlation matrix oi xiX2- In the legend, superscripts stand for the energy state and subscripts represent eigenvector contributions of 
inteipolators to make a state of interest, e.g. VpVj,Vj,v^ and V2,v_J,Vg,Vg correspond to the contributions from^i and ^2 operator for the ground state 
respectively, with larger subscript values corresponding to larger smearings. To be specific, component v" and Vj correspond to the low smearing 
sweep count of 7, v? and v" for the sweep count of 12, and so on. Errors are suppressed for clarity. 



give rise to a total spin of s - 



J, the L = 1 state 



can couple two different ways to provide a. J - ^ state. 



I.e. 



^^\ 



5 or 1 + I 



\, hence providing two 
orthogonal spin-i states in the L - 1, 70 plet represen- 
tation. Both of these negative-parity states have a width 
of «150 MeV. It is interesting that the;\fi;if2 spin flavour 
combinations successfully isolate the nearly degenerate 
ground states, in accord with the SU(6) quark model. 
This approximate degeneracy is repeated throughout 
our observed spectrum. 

In Fig. 121 similar nearly degenerate low-lying states 
from the 8x8 analysis are evident as in the 6x6 case. We 
note that, only the analysis involving;t'2 is able to extract 
a new second excited state. This state is close to the 
experimentally identified one star, N^ (2090 MeV) S \ \ 
state which has a!400 MeV of width. Further analysis is 
required with higher statistics to explore the propagation 
of these states to light quark masses. 

The eigenvectors (Fig. [8]l indicate that both xi and 



X2 operators have nontrivial contributions in extracting 
the two lowest lying states. In particular, the contribu- 
tions to the ground state is significantly dominated by 
the X2 interpolator, whereas a larger dominance of the 
X\ interpolator is evident for the first excited state, re- 
vealing a diminished role for the scalar diquark in the 
negative parity sector. This observation is in accord with 
Ref. (F 



It has already been mentioned that only the correla- 
tion matrix analysis involving ;t'2 spin-flavour combina- 
tions is able to extract the third energy state. It is evi- 
dent from the eigenvector contributions that this energy 
state is almost a pure xi state (bottom left graph). The 
fourth energy state is dominated by the x\ operator. It 
is evident that larger smearing, i.e. sweep count of 35, 
dominates over low smearing levels in isolating these 
eigen-energy states. 

The extracted masses from all the bases are very con- 
sistent. In particular, the ground state is robust. There- 
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Figure 9: (color online). Masses of the nucleon, N\ states, for all the various dimensions of correlation matrices as labelled on the upper horizontal 
axis. Dotted lines are drawn to aid in illustrating the consistency of the results. Figure corresponds to the pion mass of 797 MeV. 



fore, all the bases are considered in performing the fol 
lowing systematic error analysis. 

A systematic error analysis as that of Ref 



im 



V 



is performed to calculate the systematic errors as- 
sociated with the choice of basis, with cri, = 

N^— Tjf=i(Mi - M)2, where A^fe is the number of bases. 
This analysis is performed over all the bases consid- 
ered. The masses are averaged over the bases and er- 
rors (average statistical errors and systematic errors as- 
sociated with basis choices) are combined in quadrature. 



-4 



(^i + o-b- 

At this stage, we examine the consistency of the low- 
lying energy states over different dimensions of the cor- 
relation matrices in Fig.|9l Whereas a non-trivial mixing 
of ;^fi and xi is required to isolate the two lowest lying 
states, the near trivial mixing in the third and fourth {xi 
dominated, xi dominated) and similarly in the fifth and 
sixth states of the 8x8 analysis is manifest. 

In Fig.[TOl the negative parity ground state results are 
presented for all the correlation matrices (masses are 
given in Table [T]|. The positive parity ground state (nu- 
cleon) and Roper state results are taken from Ref. llTIl 
for comparison. Regarding the level ordering problem, 
we do not see the physical level ordering for the three 



heavier quark masses, where the A^^ state lies below 
the Roper, A^^ , state. A coincidence of the Roper and 
A^^ states occurs at pion mass of 380 MeV and the 
physical level ordering is observed at pion masses be- 
low 380 MeV. 

For evidence of the physical level ordering, we cal- 
culate the correlated jackknife error over the ratio of 
masses of Roper to the negative parity state. For our 
second lightest quark mass the result is 0.947 ^9047 and 
for the lightest quark mass it is 0.909 ^0059, providing 
strong evidence of a physical level ordering. Little ev- 
idence is found that our extracted lattice N^ state ap- 
proaches the physical state. We note that chiral physics 
is significantly altered in quenched QCD for m,, < 250 
MeV. Moreover, the finite volume of the lattice will also 
become significant in this regime. Future calculations of 
negative parity states from the variational approach ad- 
dressing the light quark masses should be done on large 
volume lattices in 2 H- 1 flavour QCD. 

At our heavier three quark masses the Nj state lies 
below the S-wave N + n state indicative of attractive 
interactions. The lower values of the N + n scattering 
state at light pion masses implies that the couplings of 
our interpolators are small for the two particle N + n 
state at these pion masses. 
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Figure 10: (color online). Mass of the nucleon, N^ , Roper N^ (Pii) and the negative parity ground state, A'j (Su) are presented. The 
combination of errors in quadrature is shown. The ground and Roper states are taken from Ref. lilVil for comparison with the negative parity ground 
state. The black filled symbols are the experimental values. 



Table 1 : Mass of the negative parity ground state of the nucleon, W ^ , for 3x3, 4x4, 6x6 and 8x8 coiTelation matrices. Masses are averaged 
over the bases described in the text and errors are a combination of average statistical errors over the bases and systematic en'ors for choice of basis. 



aM„ 


ciMll (3x3) 


aM^l (4x4) 


aMll (6x6) 


aMll (8x8) 


0.5141(19) 


1.358(19) 


1.358(20) 


1.362(21) 


1.359(23) 


0.4705(20) 


1.320(21) 


1.320(22) 


1.327(22) 


1.326(24) 


0.4134(22) 


1.279(25) 


1.278(26) 


1.285(22) 


1.283(25) 


0.3490(24) 


1.244(30) 


1.243(33) 


- 


- 


0.2776(24) 


1.223(44) 


1.221(48) 


- 


- 


0.2452(24) 


1.224(55) 


1.223(61) 


- 


- 


0.2110(27) 


1.231(40) 


- 


- 


- 


0.1905(31) 


1.249(53) 


- 


- 


- 



5. Conclusions 

An extensive analysis has been performed to ad- 
dress the issue of the level ordering problem in the nu- 
cleon sector through the variational approach. A similar 
method to that of Ref. 1 17] is used to explore this long- 
standing problem. The physical level ordering between 
the Roper and A^^ ground states in a pion mass region 
of 295 - 380 MeV is observed. However, an identifica- 



tion of a low-lying negative-parity state in the Physical 
quark mass regime remains. 
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